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Baseline calibration of a stellar interferometer is a prerequisite to data reduction of astrometric 
operations. This technique of astrometry is triangulation of star positions. Since angles are deduced 
from the baseline and delay side of these triangles, length and pointing direction (in the celestial 
sphere) of the baseline vector at the time of observation are key input data. We assume that 
calibration follows from reverse astrometry; a set of calibrator stars with well-known positions is 
observed and inaccuracies in these positions are leveled by observing many of them for a common 
best fit. 

The errors in baseline length and orientation angles drop proportional to the inverse square roots 
of the number of independent data taken, proportional to the errors in the individual snapshots of 
the delay, and proportional to the errors in the apparent positions of the calibrators. Scheduling 
becomes important if the baseline components are reconstructed from the sinusoidal delay of a single 
calibrator as a function of time. 
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I. OVERVIEW 

Astrometry as realized by a contemporary optical stel- 
lar interferometer is founded on the sensitivity to wave 
front tilts of stellar light, which leads to a dephasing pro- 
portional to the tilt and proportional to the distance be- 
tween the two telescope's apertures. This is measured 
by the amount of delay — in units of time or optical path 
length difference — added by the observatory's optical in- 
frastructure to the beam hitting the telescope closer to 
the star to adjust the phase difference of the light beams 
for close-to-coherent superposition at the detector. 

We only address geometry in this manuscript. The 
wealth of non-statistical effects of tidal Earth-crust mo- 
tion. Earth axis drifts 0, i, 1, [HI, IHIill, aberration 
and general-relativistic wavefront tilts which leave resid- 
ual noise once the known quantities are accounted for is 
left aside. 

Mathematics of small differences demonstrates in Sec- 
tion |TT] how measurement of the difference in the delay 
defined by an angular separation of two "science" ob- 
jects in the celestial sphere requests some knowledge of 
the baseline vector — which we vaguely define as the sepa- 
ration between the two input pupils of the two telescopes 
involved. Baseline calibration is the auxiliary observation 
of well-known (in the astrometric sense) calibrator stars 
to deduce this baseline geometry. 

On can think of two pure forms. First there 
is a Fourier mode which matches the sinusoidal de- 
lay as a function of time with the three free pa- 
rameters of the baseline vector tracking an individual 
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star (Section IIIip . This minimizes the time overhead 
of slewing/pointing/acquisition/fringe-loop-lock cycles. 
A characteristic sensitivity of the baseline parameters 
means that is is favorable to reserve time slots six hours 
apart [18]. 

This compares to the second form — which one may call 
field mode — in which delays of a larger set of calibrator 
stars covering a wide range of altitudes and azimuths are 
gathered (Section llV|) . 

The topic of this script is focused on the narrow ques- 
tion: given statistical errors of the delay measurement in 
Fourier mode and statistical errors in positions derived 
from star catalogs in the field mode, how long or how 
many of them, respectively, does the calibration need to 
balance their effect to the levels set by the astrometric 
mode. 



II. FUNDAMENTAL TRIGONOMETRY 
A. Daily OPD 

We start with a summary of the fundamental geom- 
etry of observing a star with perfectly stable telescopes 
orbiting a fixed Earth axis without atmospheric or simi- 
lar distortions [lj|. 

The terrestrial coordinate system defines geographic 
longitude A, latitude (p and altitude H of two telescopes. 
By conversion into a Cartesian frame and construction 
of the mid-point between any pair of these, we define 
the geographic longitude A and geographic latitude (f) of 
the baseline, which serves to define the topocentric Alt- 
Az system for the baseline — which implies that af, = 
since we define an individual system for each pair. In 
this topocentric coordinate system, the baseline vector 
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TABLE I: A model of telescope coordinates and the associated 
equatorial variables of the Very Large Telescope Interferome- 
ter aligned with WGS84 conventions. 



telescope 




A (rad) 


(f) (rad) 


H (m) 


Ul 




1.228800386 


-0.429833092 


2635.43 


U2 




1.228796107 


-0.429825122 


2635.43 


U3 




1.228790929 


-0.429819523 


2635.43 


U4 




1.228780856 


-0.429823005 


2635.43 


baseline 


b 


A 




Sb ht 




(m) 


(deg) 


(deg) (deg) (deg) (deg) 


U12 


56.4 


-70.4050 - 


24.8342 -153.9 


46.8 -49.7 


U13 


102.1 


-70.4048 - 


24.6272 -147.5 


43.9 -56.8 


U14 


130.2 


-70.4045 - 


24.6273 -119.5 


25.6 -76.8 


U23 


46.5 


-70.4047 - 


24.6270 -139.8 


39.8 -63.8 


U24 


89.5 


-70.4044 - 


24.6271 -98.6 


7.8 -86.4 


U34 


62.5 


-70.4042 - 


24.6269 -69.3 - 


-18.4 -99.0 



1 sid. day, 2 tz 




FIG. 1: Equation ^ generates a periodic delay in time. 



has length b and splits into Cartesian coordinates 

- cos Ab cosflb 
sin Ab cos flb 
sin ah 



(1) 



as a function of azimuth A}, (South over West) and incli- 
nation af,. A rotation matrix U 



U = 



— sin (f> cos A 

sin A 
cos 6 cos A 



— sin (f> sin A cos < 
— cos A 
cos (h sin A sin ( 



(2) 



transforms these coordinates to Cartesian coordinates in 
a geocentric frame for suitably defined celestial declina- 
tion Sb and hour angle ht p^ . 



h = U -b 



cos Sf, cos(A — rih ) 
cos Sb sin(A — hb) 
sin 6b 



(3) 



Table H] illustrates this transformation for a Chilean site. 
(The constant height H above the ellipsoid implies that 
neither telescopes nor baselines are co-planar.) 

The direction to the star of right ascension a and dec- 
lination S at current azimuth A and zenith angle z is 

-cos^sinzN / cos (5 cos(A — /i) 

sin A sin z = ■ cos (5 sin(A — /i) | . (4) 
cosz J \ sin (5 

The equation-of-motion of the optical path difference 
(OPD) £) is [13 

D = s h = 6 [cos z sina;, -I- sin z cos flb cos(^ — ^b)] 

= b [sin S sin Sb + cos S cos Sb cos(/i — hb)] .(5) 

This is the common spherical coordinates formula for the 
angular distance between two points, one fixed at (a, S), 



the other at fixed Sb cycling the polar axis in one sidereal 
day Td, at an angular velocity of 



LOd — 2TT/Td « 73/irad/s. 



(6) 



The characteristic parameters of D{h) are sketched in 
Figure [1] a time-independent offset 5 sin 5 sin <5h , an am- 
plitude b cos S cos Sb, and a phase hb. The baseline length 
b and the angles Sb and hb are encoded in phase and 
amplitude of D plotted over time t. 

The time-dependent projected baseline angle pb is the 
position angle of high sensitivity of D to changes in the 
sky coordinates, that is the direction of high interfero- 
metric resolution UM. 



t&n Pb = 



cos Sb siii{h — hb) 



cos S sin Sb — sin 6 cos Sb cos(/i — hb) 



(7) 



B. Differential OPD 

We are concerned with differential astrometry of ob- 
serving two objects at a single point in time rather than 
switching between the two stars P,[2ll. We note ^ for 
two stars with coordinates Si and ai = LST — hi , mean 
positions S and a, 

Aa = a2-ai; a = (q;2 + ai)/2; (8) 
= S2-Si- S={S2^Si)l2. (9) 
h = {h2 + hi)/2 = LST - a. (10) 

The time-independent cosine of the angular separation t 
is 

COST = cos(5i cos(52 cos(q;i — 012) + sin(5i sin^2 (H) 

1 



l-i(A5)2-icos2 5(Aa)2 



24 



{\sY 



+ -(A5)2(Aa) 



1 

24 



cos^ ^(Aa)-* -t- • ■ • (12) 
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FIG. 2: Equation (|13|l defines a differential delay with the 
same period length as the wide-angle delays Di and D2- 



in these coordinates, neglecting sixth and higher order 
mixed differentials. The differential OPD becomes a si- 
nusoidal function of time, too, 

AD = D2-D1 = (s2-si)-b = o+pcos(h-hb+ho), (13) 

which defines a coordinate offset o, a daily amplitude p, 
and a time shift ho (Fig. [2]). For small Aa and small AS, 
up to fourth order in A6, 



o « bAS cos S sin ( 



1 

24 



b{A6f cos 6 sin 6b. (14) 



Up to fifth mixed order in the differentials we have the 
squared amplitude of the differential delay. 



p2 « cos^ 6b[sm^ 6{ASf + cos^ d{Aa) 
-^cos'SiAaf-^iASfiAaf], 
and the shift in hour angle 



1 

12 



sin^ 6{AS)^ 
(15) 



tan hn « cot S—^ [l 
AS 



12 



(AS)' 



12 



(Aa)' 



1 

'720 



(Asy 



1 

144 



(A<5)2(Aa)2 + — (Aa)^]. (16) 



The right hand side contains a factor cosSAa/AS, the 
tangent of the position angle of the binary to lowest order 
in the differentials. 



the sky, which could either be represented as Aa and AS, 
or alternatively as distance r and position angle. 

From a fit to a sine model like a measurement can 
extract three parameters, which represent the six free 
Cartesian components of head and tail of the baseline 
vector minus the degrees of freedom of the three compo- 
nents of the baseline center (which represent a free rigid 
translation of the baseline and do not contribute to the 
interferometric signal) . 

For the purpose of this script, the time base is not con- 
sidered an independent source of error. (On a real-time 
bus an individual query for a time stamp has a resolu- 
tion of roughly 10 fxs, i.e., 150 fias after multiplication by 
LUd- Micro-controller boards usually govern the readout 
process, so the jitter is much smaller.) 

These three parameters could either be stored as 
Cartesian coordinates of the baseline vector, or in the 
coordinates 6, S^ and hb which are more meaningful in 
the context of the analysis of D{t) in general. 

Independent of this question of format, there is some 
redundancy between the three fitting parameters o, p and 
ho contained in D{t) and the two positional parameters 
of the binary, supposed an independent baseline calibra- 
tion provides the auxiliary b, 6b and hb- The task of the 
baseline calibration is to support inversion of equations 
P^ - pB|) . From this point of view, we need the product 
b sin Sb to reduce pi)) . the product b cos Sb to reduce (fT5|) . 
and hb to reduce ([16]). 

If one of the three equations is not activated for some 
reason, one of these three projections of the baseline vec- 
tor (onto the polar axis and on the equatorial plane) does 
not need to be calibrated either, because — in principle — 
two equations for two unknowns remain. The U24 base- 
line in Table [J provides an example of this aspect, where 
b sin Sb ~ 12 m. The factor cosi5 might reduce the prod- 
uct b cos 6 sin 6b on the right hand side of to only 6 
m, which couples an accuracy of 100 ^as= 5 x 10^^" rad 
in AS to an accuracy in o of 3 nm. If the astrometric 
run cannot meet this requirement, equation p4|) drops 
out of the data reduction process, and in turn, the value 
6 sin (5b is not requested from the baseline calibration. If 
the value of b cos 6b is very small, equation (|15p may be- 
come disposable with the same rationale. 

From a similar mathematical but entirely different en- 
gineering point of view, the measurement of D will likely 
be assisted by a metrology system which is difficult to cal- 
ibrate over longer periods of time (which includes drifts 
of air densities if operating in air and atmospheric lensing 
effects) and is difficult to run without interruption. The 
time derivative, i.e, the delay velocities 



Correlation of Variables 



AD{t) ^ -pujdsin{h - hb + ho), (17) 
D(t) = —bu!d cos 6 cos Sb sin(/i — hb) (18) 



The astronomer's interest of astrometric operations is will be available more easily, but this eliminates the off- 
in the measurement in the two parameters that charac- set, including the differential offset o and (jl4p . as an in- 
terize the relative position of the binary's components on dependent piece of information. 
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D. Requested Baseline Accuracy 

The "cushion" and "barrel" distortions of the geometry 
by the cubic terms of (fT4)) and biquadradic terms of (fT5)) 
and (fT6|) are usually negligible: If the angular distance is 
limited by some finite field-of-view to t < 2', the squares 
are limited to < 3 x 10"'' rad^. If the other terms 
in the sum are of the order of unity, these higher order 
contributions change A6 or Aa by less than 2' x 3 x 
10-^ « 40 //as. 

So requirements on the baseline calibration can be 
derived from the first orders o ~ bAS cos S sin Sb and 
p ~ 6 cos (Sb sin (5 A(5 ~ 6cos(56COS(5Aq;. In a broad sense 
the relative error in AS and Aa is of the order of the rel- 
ative error in o plus the relative error in b sin 6b, or p plus 
the relative error in b cos Sb, respectively. As a guideline: 
an accuracy of 100 /las in a field of t < 2', or 10 /xas 
in a field of r < 10", is equivalent to a relative accu- 
racy of 1 X 10^®, which splits evenly into requirements of 
5 X 10^^ in o (p) and b, which is a requirement of 50 /im 
for a baseline of 6 = 100 m of length. 

Equivalent considerations are available for p^ : The 
relative error in the value of b matches the relative error 
in the value of D. For an error 50 /^m in a 100 m baseline 
we need a relative accuracy of Z) of 5 x 10^^. Since the 
speed is of the order 1 cm/s for a 100 m baseline, we 
want to have D in absolute precision to 5 x 10~^ m/s = 5 
nm/s. Near the turning points h « /if,, the requirements 
for absolute precision become tighter, enforced by the 
sine factor in (fT5|) . 

The signal D{t) may also be small because either S or 
6b are near 90°, which puts higher stress on the measure- 
ment of the absolute value of this velocity. The effect 
through cos 6 means that a star on a closed apparent 
orbit around the corresponding celestial pole of the ob- 
servatory's hemisphere induces a small variation of de- 
lay because its apparent position does not change much 
anyway. The effect through cos 6b punishes North-South 
baselines of observatories near the equator because 6b 
measures the angle between the baseline and the equato- 
rial plane. 

To keep the number of parameters basically down to 
three, this manuscript does not introduce models of tele- 
scope axes runouts as a function of azimuth and altitude. 
They do have an impact on the calibration: a runout of 
50 fim by one of the telescopes, i.e., not matched by the 
other, is equivalent to that characteristic baseline tilt of 
100 mas over b = 100 m mentioned above. 



III. TRACKING A SINGLE STAR 

A. 3-parametric fit 

The simplest scheme to derive baseline coordinates 
from reverse wide-angle astrometry is to monitor the 
trace for a single star, as in Fig. [2 over a significant 



portion of one of the periods. 

One reason for such a strategy can be that there are 
few calibrator stars with small errors in their (a, 6) coor- 
dinates and mixture with data from more stars — strategy 
of Section IIVI — could deteriorate the statistics. Another 
reason is enhanced efficiency if the calibrator star is one 
of the components of the binary system; baseline calibra- 
tion and astrometric observation may then run concur- 
rently if delay and differential delay data are recorded 
simultaneously. 

A Fourier analysis can extract amplitude, offset and 
phase by a least-squares fit 

N 

[b sin 6 sin Sb + b cos S cos 6b cos{hj — hb) — Dj]'^ min 

(19) 

by N measurements. The individual delay measurements 
Dj are characterized by an error variance a^{D). One 
might include some statistical weights of the squares 
in (|19p by some inverse function of the projected 
baseline — which represents the angular resolution of the 
interferometer — but this adds more parameters to the 
present analysis and obscures the principles laid out be- 
low. 

The output (and free parameters) of the fitting pro- 
cess are b sin 6b (polar component of the baseline) , b cos 6b 
(equatorial component of the baseline) and hb (interfer- 
ometric hour angle). The error propagation of such a 
Fourier analysis has been discussed in the literature [H, 
§7.04] [l^. There is no need to split the components 
into length b and angle Sb for two reasons: 

1. These are correlated variables. There is no profit 
from dividing the information (plus tracking the 
covariances), to re- unite them afterwards for appli- 
cation in 111]) and P^ . 

2. If the the observable is the delay velocity D and the 
differential delay velocity AD as discussed above, 
the components remain properly separated. One 
can basically ignore the error analysis of the offsets 
(polar components). 

B. Scheduling 

The limits of observing during the night and observing 
stars at some minimum altitude above the horizon lead 
to a characteristic error propagation depending on which 
portion of the delay is covered by the observation. As 
sketched in Fig. [31 N observations of snapshots Dj would 
start at some time ifo after the daily maximum. The 
quality of the fit to the amplitude or offset depends on 
whether ip is near zero or tt on one hand, or 7r/2 away 
from a maximum or minimum on the other hand. 

We summarize simple Monte Carlo calculations of 
these dependencies in Figures |3HS] for a baseline of 6 = 
100 m. The common feature of Fig. 2] and O is: 
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FIG. 3: The baseline calibration in Fourier mode samples 
the delay of Fig. [1] at points j = 1, 2, ... A*' in time, starting 
at ipo- 
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FIG. 4: Errors to the delay oflFset, delay amplitude and de- 
lay hour angle for A*' — 360 measurements in intervals of 10 
seconds or N — 180 measurements in intervals of 20 seconds 
(an observation taking one hour). (j{D) — 1 fj,m for each 
measurement. 



• The quality of the measure of the offset (mean) does 
not depend on when the observation is started. 

• The amplitude of the curve is obtained with highest 
quahty if the observation covers the part of maxi- 
mum delay velocity. To achieve highest precision in 
the hour angle, however, the observation ought to 
cover the part of zero velocity near an extremum, 
six or eighteen hours earlier or later. 
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FIG. 5: Errors to the three parameters of Fig.|3]for A'' — 180 
measurements in intervals of 10 seconds, covering only the 
next half an hour after ipQ. 
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FIG. 6: Errors to the three parameters for measurements at 
intervals of 10 seconds, total integration time from 0.1 to 1 
hour, a{D) = 1 nm or 0.5 fim, (po = 2 hrs. 



Fig.[4]shows the generic dependency of errors ^ N^^^"^ if 
duration and time slot of the observation stay the same. 
The transition from Fig. |4] to Fig. [5] demonstrates that 
cutting the observation time by half, keeping the detec- 
tor integration time the same, increases the errors by an 
approximate factor of six. Fig. [S] emphasizes this depen- 
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0.4 as, 100 nm 
0.2 as, 100 nm 
0.1 as, 100 nm 



0.05 as, 100 nm 
0.025 as, lOOnm 
0.025 as, 2 (im 



Figure [8] uses a more constrained region of the sky 
with z < 40° and achieves inferior accuracy for equivalent 
numbers of stars. (The side effect that the projected 
basehne is larger on the average than in Fig. [7] which 
implies better resolutions is not taken into account.) 



100 




10 



FIG. 7: The error a{b) to a baseline length b = 100 m from 
a least-squares fit to A'^ measured delays D for A'' stars. 



dence of the errors '--^ A^^5/2 ^^^^ constant detector inte- 
gration time), but confirms the expected proportionality 
to the error a{D) of the individual readout. 

In this model, the error <t{D) is an effective superposi- 
tion of an error induced by the error in the interferometric 
phase plus an error from the jitter in the time base. At 
velocities D < 1 cm/s, an error of cr(ft.) = 20 /iS in the 
time base is equivalent to an error (7{D) = Da{h) < 200 
nm, for example. Techniques to reduce this error by im- 
plementing detector-readout schedules in low-level micro- 
controller programs are not in the scope of this paper. 



IV. CALIBRATOR STAR CATALOG 
IMPRECISION 

A. Baseline Length Calibration 
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FIG. 8: The errors of Figure [7] if baseline calibrator stars are 
only selected in the zenith range z < 40°. 

The apparent positions of stars are effected for exam- 
ple by the chromatic dispersion in air, the product of 
the dielectric susceptibility of air [l^ at the telescope — 
measured in radians — and the tangent of the zenith an- 
gle. Examples for two infrared windows for Paranal con- 
ditions are plotted in Fig. [9] and [10] for three different 
water vapor densities. 

This implicit spread of a few tens of mas depending on 
the star color is of the same order as extrapolation errors 
for the apparent position from integrated proper motions 
since the Hipparcos epoch [l3| ■ 



The estimated error a{b) in the baseline length for b — 
100 m obtained by a 1-parametric least squares fit after 
visiting N stars on the sky that are rather homogeneously 
distributed in the range z < 60° is shown in Figure [T] 
The N positions are randomly selected over the sky in 
the zenith range z < 60°, taking subsets of the Hardin- 
Sloan-Smith points [?!, and have been displaced by angles 
with five different Gaussian widths between 0.025" and 
0.4" to generate a measured delay, to which in addition 
errors of a{D) = 100 nm or a{D) = 2 /um are added. 
Double logarithmic axes scales are chosen to verify that 
the reduction in a{b) is approximately proportional to 

iV-l/2. 



B. Baseline Vector Calibration 

Fig.[TT]are results of a least squares fit of the 3 degrees 
of freedom of the baseline vector to a series of delay mea- 
surements, minimizing ^i[Dj — Sj ■ b]^ over a set of 
"noisy" star positions Aj and Zj and building a statis- 
tics over the variables 6, Af, and ab- The error to the 
the measured delay has been set to zero, because Figs. 
[THI] reveal that they are not important if they remain 
a{D) < 2 fim. The N positions have been rather ho- 
mogeneously distributed over the sky as in Fig. [7] in the 
zenith range z < 60°, and have been randomly displaced 
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FIG. 9: The chromatic difFerential transverse atmospheric dis- 
persion (TAD) in the H band for tanz = 1 at an ambient 
pressure of 744 hPa, relative to a wavelength of 1.65 fim. 
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FIG. 10: The differential transverse atmospheric dispersion 
in the K band for tanz = 1 at an ambient pressure of 744 
hPa, relative to a wavelength of 2.2 ^m. 



by angles with different Gaussian widths of 0.4" (pluses), 
0.2" (crosses) or 0.1" (stars). 

The accuracy in the baseline length b is the same as 
obtained with the 1-parameter fits of Fig. [T] Clearly, for 
small N, the error in the angles just echoes the error 
introduced in the star positions. 

The elevation angle at (tilt of the baseline versus the 
horizontal) is obtained roughly twice as accurate as the 
azimuth Af,. The interpretation of this bias is: the even 
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FIG. 11: The error to a baseline length b = 100 m, and to the 
baseline orientation angles from a 3-parametric least-squares 
fit measured to delays D of N calibrator stars. 



distribution of the calibrator stars along azimuths and 
their more clumpy, overhead distribution along zenith 
angles means that a measurement via the projections D 
achieves low resolution along the horizontal for the two 
subsets of stars in the two opposite pointing directions of 
the baseline (TbI, Fig 2-8] . The coupling along the vertical 
coordinate is stiffer on the average, and the information 
contained in the delays better distributed to deduce the 
baseline tilt (pitch) versus the horizon than the baseline 
rotation (roll) around the zenith. 

This difference in fitting quality in the {ab,Af)) angles 
is more obscure in the {6b, hj,) system, because multipli- 
cation with the inverse (transpose) of the JJ-matrix ^ 
weights components depending on the geographic lati- 
tude (f). Anyway, the angles in Fig. [11] are related to 
the effect of Earth axis pointing [8| , which can be writ- 
ten down as a change of the effective {Sf,, hi,) coordinates. 
Formula ([5]) is in fact symmetric as one could swap the 
variables S and 6b or hj and hb without changing the 
delay. Clearly, the calibration measures a baseline orien- 
tation in an ecliptic (celestial) reference frame, not in a 
terrestrial reference frame. The two sides of this coin are 

• Pushing errors in af, or Ab below the lim- 
its set by the errors in the lERS angles 
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www . usno . navy . mil/USND/earth-or ientation 
/eo-products is useless if the aim is to reduce 
delays to the terrestrial baseline data. 

• There is a prospect of setting up a competitive opti- 
cal reference frame to 40 mas — aside from the influ- 
ence of all systematic effects — if iV > 15 stars with 
positions accurate to 100 mas are available. (This 
estimate follows almost trivially from a N~^^^ scal- 
ing of uncorrelated errors.) 

A further remark: The baseline vector calibration is in 
a general sense equivalent to determining the geographic 
longitude and latitude of either head or tail of the base- 
line vector: one can change the orientation of the baseline 
by either tilting it explicitly or keeping it always hori- 
zontal and sliding it with the tangent plane across the 
Earth surface. In mathematical prose: Introducing the 
geographic latitude (j) or longitude A as additional fitting 
parameters into the minimization procedure defines an 
ill-conditioned problem. 

V. SUMMARY 

Baseline calibration reduces measurements of the 
scalar variable of the delay to baseline vector compo- 
nents along the polar and equatorial axes, and splits the 
equatorial component into two with the aid of clocks. 

The interest is in the measure of angles, such that the 
requirements on lengths (delay and baseline) are not for- 
mulated in absolute but in relative units: the relative 
error in the star separation of the astrometric observa- 
tion is the relative error in the delay measure plus the 
relative error in the baseline length, and similar generic 
statements apply separately to the components in right 
ascension and declination. 

The statistical errors in the baseline coordinates are 
proportional to the errors of the individual delay (or its 
speed) and time stamp, and proportional to the inverse 
square root of the number of independent measurements, 
as expected. If the baseline coordinates are derived from 
a Fourier analysis of the delay or delay speed as a daily 
function of time, the statistical errors depend more deci- 
sively on schedules and on the percentage of the sidereal 
period that is covered. 
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APPENDIX A: TIDAL MOTION 

One contribution of the definition of the telescope coor- 
dinates in an extra-terrestrial coordinate system is given 



by the influence by the ocean tides that load and release 
the non-rigid Earth crust According to the GOtOO . 2 
model by Bos and Scherneck 0, the amplitude for the 
Paranal geographical coordinates is < 2 cm in vertical 
and < 0.7 cm in horizontal directions: Fig. [121 
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FIG. 12: Illustration of the vertical and lateral motion of 
the Earth crust at the Paranal coordinates parametrized with 
11 harmonic frequencies 0]. The fundamental mode exhibits 
the two daily tides, the envelope is governed by the Moon's 
period. 



The Earth crust motion in other regions of the planet 
may be up to 10 cm. 

The periodic influence on astrometry can be estimated 
by converting the East- West motion into a change of the 
instantaneous geographic longitude A, the North-South 
motion into a change of the geographic latitude (f). A 
sliding by 1 cm translates into a tilt of 0.01/6.38 x 10^ ~ 
1.5 X 10"^ rad w 0.3 mas, and is therefore not relevant 
for the baseline calibration if the errors in calibrator star 
positions are roughly one or two magnitudes larger. With 
a similar rationale, a change of 5 in (fT4|) or (fTS]) by 1.5 x 
10~^ remains negligible if the request for relative errors 
in o or p is of the order 5 x 10~^ as argued in section Hi D I 
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